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Abstract. The three-dimensional inviscid Hasegawa-Mima model is one of 
the fundamental models that describe plasma turbulence. The model also ap- 
pears as a simplified reduced Rayleigh-Benard convection model. The math- 
ematical analysis the Hasegawa-Mima equation is challenging due to the ab- 
sence of any smoothing viscous terms, as well as to the presence of an ana- 
logue of the vortex stretching terms. In this paper, we introduce and study a 
model which is inspired by the inviscid Hasegawa-Mima model, which we call 
a pseudo-Hasegawa-Mima model. The introduced model is easier to investi- 
gate analytically than the original inviscid Hasegawa-Mima model, as it has a 
nicer mathematical structure. The resemblance between this model and the 
Euler equations of inviscid incompressible fluids inspired us to adapt the tech- 
niques and ideas introduced for the two-dimensional and the three-dimensional 
Euler equations to prove the global existence and uniqueness of solutions for 
our model. Moreover, we prove the continuous dependence on initial data of 
solutions for the pseudo-Hasegawa-Mima model. These are the first results 
on existence and uniqueness of solutions for a model that is related to the 
three-dimensional inviscid Hasegawa-Mima equations. 
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1. Introduction 

The three-dimensional inviscid Hasegawa-Mima equations were first derived by 
Hasegawa and Mima as a simple model that describes plasma turbulence. The 
Hasegawa-Mima equations are also one of the simplest and most fundamental mod- 
els that describe the electrostatic drift waves. The three-dimensional equations that 
describe the coupling of the drift modes to ion-acoustic waves that propagate along 
the magnetic field are given by the system (c.f., [1], [15], [18], [38]): 

^(A^0 - 0) + J(0, A„0) + - ^ = 0, (1.1a) 
at oy oz 

^ + J(0,..) + |^=O, (1.1b) 
where J is the two-dimensional Jacobian defined by J(f,q) — tt-tt ~ tttt and 

J \j 1 ^ ; ay oy ox 

A/i = + is the horizontal Laplacian. Here (j) is the electrostatic potential 
and is simultaneously the stream function for the horizontal flow in the xy-plane, 
so that the horizontal velocity is given by = z x Vh(/>. Here Vd is a constant 
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that is proportional to the density gradient and Vz is the normahzed (by the ion- 
acoustic wave speed) velocity in the z-direction. Notice that the last term in the 
first equation is the "ghost" of the vortex stretching. 

In the context of geophysical fluid dynamics, the three-dimensional inviscid 
Hasegawa-Mima equations appear as a reduced Rayleigh-Benard convection model 
which describes the flow motion of a fluid heated from below. The three-dimensional 
rotationally constrained convection model was derived in [21], and [29], under 
the assumption of very small Rossby number. In the three-dimensional domain 
Q, ~ [0,X]'^, the equations are: 

^ + (u . V.).. - Mz^^OlM ^ ^ i ^^^^ ^^ ^a) 

at oz Re 

^^ + (u.V,0^-- = -A.., (1.2b) 

f: + (u.v.)e' + .f = i^A.e', (1.2c) 



1 a^e d[Q'w) 



= 0, (1.2d) 



Pe dz"^ dz 

e' dxdy = 0, (1.2e) 

[0,L]2 

Vft-u = 0, (1.2f) 

where the horizontal-average of a given function is defined as: 

1 



L2 



(j) dxdy. 

[0,L]2 



Here, u = (u,w) is the horizontal component of the velocity vector field {u,v,w), 
and the vertical component of the vorticity, uj ~ V/i x u := ^ — and the 

fluctuation of the temperature, Q = — G), are unknown functions of {t;x,y,z), 
while the horizontal-mean temperature & is an unknown function of (t;z). Re 
is the Reynolds number, F is the buoyancy number, Pe is the Peclet number, 
Vh = (^4)and A, = ^-K|^. ^ 

The two-dimensional Hasegawa-Mima equations were derived by Charney and 
Obukhov as a shallow water model from the Euler equations with free surface under 
quasi-geostrophic balance assumption. The geophysical model that was obtained 
is: 

^^^"'^^r + Ji^o,^hM + J(</>o>i3 + M = 0, (1.3) 

where 4>o{x, y) is the amplitude of the surface perturbation at the lowest order in the 
Rossby number, (I>b{x) is a given function that parametrizes the bottom topography, 
and here again J is the Jacobian while F is the Froude number constant. When 
a flat bottom is considered; that is when 0b is taken to be equal to constant one 
obtains 

\~ J{(P0,^h(P0)+ P-g^ ^0: (1-4) 

which is the inviscid Hascgawa-Mima-Charney-Obukhov (HMCO) equation, where 
00 = 0o(^;3:, y). Global existence and uniqueness of solutions of the inviscid 
HMCO equation can be easily established following the same ideas as for the two- 
dimensional Euler equations. Indeed, it was remarked in [10], and in [26], that (1.4), 
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with F ^ /3 ~ 0, is equivalent to the two-dimensional Euler equations of incom- 
pressible inviscid fluid, where (jjo plays the role of the stream function, w = A^^o 
the vorticity and u = (— is the velocity field. The existence and unique- 
ness of strong local solutions of (1.4) in iJ*(K^), with initial data (/fp G 
for s > 4; and the existence of a weak global solution in i?^(E^), with initial data 
00 S were established in [26]. The uniqueness of a global strong solution 

in 7J*(K^), for s > 4, was later established in [10]. In a recent paper [8], we followed 
the work of Yudovich [34] to prove the existence and uniqueness for the vorticity 
A/100 in the space in the case of no- normal flow boundary conditions, i.e. 
Dirichlet boundary conditions for the stream function 0o • Notice that 0o G H'' , for 
s > 4. is equivalent to Ah4>o E H'^~^. Thus, our result, reported in [8], improves 
the previous results, since Ah(po G H''~^ C L°° for s > 3. 

The two-dimensional viscous version of the Hasegawa-Mima equations was also 
studied, later on, by some authors using standard tools used for studying the two- 
dimensional Navier-Stokes equations. For example, the existence of a global at- 
tractor in H'^{n) of the viscous two-dimensional equations, was studied in [36] and 
in [37]; and an upper bounds for its Hausdorff and fractal dimensions were also 
established. A perturbed version of the two-dimensional Hasegawa-Mima equa- 
tions was studied in [11], the authors considered the two-dimensional Hasegawa- 
Mima equations with a hyper-viscosity term and used techniques introduced for the 
Kuramoto-Sivashinsky equation (KS) by Nicolaenko, Scheurer and Temam [25], 
for odd-periodic solutions. The authors in [38] studied the existence of a global 
attractor of a viscous three-dimensional model similar to the three-dimensional 
Hasegawa-Mima equations that they called a "generalized" version of the three- 
dimensional Hasegawa-Mima equations and they established upper bounds for its 
dimension. We emphasize that here the authors followed well-established tools for 
dissipative equations (see, e.g., [32] and reference therein) to prove the existence 
of the global attractor and estimating its dimension. The added artificial three- 
dimensional viscosity simplifies the mathematical analysis drastically and allows 
for the implementation of these well-established tools. We stress here that all the 
previous studies have been focused on the viscous and the inviscid two-dimensional 
Hasegawa-Mima equations, and on the viscous three-dimensional case following 
well-established tools. 

Investigating the global existence of (1.2) is challenging due to the fact that the 
unknown functions depend on the three spatial variables (x, y, z), while the regu- 
larizing viscosity terms acts only in the horizontal variables {x,y). The situation 
for (1.1) is even more challenging since it has no viscous regularization at all. In- 
spired by these models, i.e. (1.1) and (1.2), we introduce the following "simplified" 
mathematical inviscid model - the pseudo-Hasegawa-Mima equations: 

^ + (u • V,)«; - (7oif^ - 0, (1.5a) 
at oz 

- + (u.V.).---=0, (1.5b) 

Vh-u = 0, (1.5c) 

for some Uq velocity scale in the z-direction. Here u — (u, v) is the horizontal 
component of the velocity vector field {u,v,w), and the vorticity uj = x u := 
^ — 1^. One can see easily that (1.5) is a simplified version of the Hasigawa-Mima 
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equations (l-l)- The goal of introducing, and investigating, this mathematical 
model is to shed hght on the analysis of (1.1) and (1.2). The system (1.5) is 
simpler because it has a nicer mathematical structure. Specifically, let us denote: 

O — w + Lui, rj — w^Luj. (1-6) 

Adding and subtracting (1.5a) and (1.5b) yield the coupled system: 

^ + (u.V,)0-C/o|^ = O, § + (u.V,0r? + (7o^ = 0, V;,.u = 0, (1.7) 
ot oz ot az 

where V/i x u = w = ^{9 — rj). System (1.7) is similar in its structure to the 

two-dimensional Euler equations due to the relation between u and the vorticity w, 

where the variable z is thought of as a parameter, and also in the sense that it is a 

purely three-dimensional transport system that does not have a vorticity stretching 

term. Notice that the velocity fields that transport 9 and rj in (1.7) are not the 

same. This is in fact a rcfiection of the vortex stretching terms in the original model 

(1.5). 

We denote by T'' the L-period box [0,i]'^. We will establish in this work the 
global existence and uniqueness for system (1.7) in T'^, i.e. subject to periodic 
boundary conditions, and then use this result to conclude the global existence 
and uniqueness of solutions for system (1.5) in T'^. As in the case of the two- 
dimensional Euler equations in vorticity formulation (see, e.g., [23] and [24]), we 
have the analogue of the two-dimensional periodic Biot-Savart law: 

u{x,y, z) ~ K * io -.^ / K{x — s,y — ^) io{s,^, z) dsd(,, (1.8a) 

h Jj2 

K{x,y)^WiG{x,y), (1.8b) 

where G(a;, y) is the fundamental solution of the Poisson equation in two-dimensions 
subject to periodic boundary conditions, the binary operation * denotes the hori- 

h 

zontal convolution, and = (l-^) '^^ explicitly restrict ourselves 

to the unique solution u of the elliptic system: V/i x u = to and V/i • u = 0, that 
satisfies the side condition u(2;, J/i z) dxdy = for every z e [0, L]. 

In section 2, we state and prove a technical two-dimensional logarithmic inequal- 
ity that improves the Brezis-Gallouet [2] and the Brezis-Wainger [3] logarithmic 
Sobolev inequalities. This logarithmic inequality allows us to prove some useful 
estimates that we use in proving the global existence and uniqueness results. In 
section 3, we will study the global existence and uniqueness of smooth solutions of 
a "Hnearized" version of (1.7) and establish some relevant uniform estimates. In 
section 4, we use the established results and the uniform estimates obtained for the 
"linearized" system, in section 3, to prove a global existence and uniqueness result 
of solutions, and obtain some uniform estimates for (1.7). Even though system 
(1.7) does not seem to have explicit vorticity stretching term, but the terms ^ and 
^ can be thought of as if they are vorticity stretching terms and produce similar 
challenge. The main challenge in this work is in obtaining the relevant uniform 
bounds that we proved based on estimates we obtain in section 2, using the log- 
arithmic Sobolev-type inequality. We stress that our results are one of the first 
results related to the inviscid three-dimensional Hasegawa-Mima equations. 
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2. Preliminaries 

In this section, we introduce some preliminary material and notations which are 
commonly used in the mathematical study of fluids, in particular in the study of 
the Navier-Stokes equations (NSE) and the Euler equations. 

Let T be the set of all trigonometric polynomials with periodic domain T'^ := 
[0,i]3. We denote by iJ'(T3), ^/^^^(T^), H^iJ^) = W'^iT^) to be the closures of J" 
in the usual Lebesgue and Sobolev spaces, respectively. Since we restrict ourselves 
to finding solutions over the three-dimensional L-periodic box T'^, therefore, we 
can work in the spaces defined above consistently. We define the inner product on 
L2(T3): 



(u, v) = / u(x) • v(x) dx, 

JT3 



|^2^,j,3-, = (u, u)^^^. (We use this notation for both 
scalar and vector fields). We define the following Banach spaces: 



LP(T3) := |/ : ||/(a.,y)!|^5(T3) := \f{x,y.z)\'' dz < oo, for every (x,2;) S [0,L]'j 
LliT') := |/ : 11/(^)11^^(^3) := £ £ |/(a^,2/,^)r dxdy < oo, for evcryz S [0,L]| 
L'f^iT^) < / : I1/(^:)I1l=c(t3) := esssup \f{x,y,x)\ < oo, for everyz e [0, L] \ , 

LTfi^^) ■= |/ : II/IIl»(l^(T3)) ■=^|^j^"P II/(^)IIl^(T3) < ooj ■ 

Similarly, we define the spaces L^{W^'P{T^)) and Lf{H^{T^). We also denote by 
the Banach space: 

l.P{T^) := LP(T3) X L'PiT^) 

and its associated norm: 

II('A»IIlp(t3) = !1'?!'I1lp(t3) + I1V'I1lp(t3) • 

Similarly, we define the spaces Lf (T'^), LP(T'3), L^(T'3) and L^^''(T3). The natural 
definitions of |l-|li4/-'"'i'(T3) ^^'^ \\-\\w'"-''{t^) ^'-'^ 

II, ,11^ . \ ^ 7- ' „ I l7-ia../_ ..M^ 

(T^) • — 

|a|<m 

and 

|1u(z)|12^,„,p,t2) :- V / \D^n{x,y,zf dxdy. (2.2) 

\a\ <m 

Let F be a Banach space. We denote by LP{[0,T];Y) the space of (Bochner) 
measurable functions t i-)- w{t), where w(t) € Y, for a.e. t € [0,T], such that 
the integral J^^ \\'w{t)\\Y dt is finite. A similar convention is used in the notation 
C'^([0,T];F) for fc-times differentiable functions of time on the interval [0,T] with 
values in Y. We denote by C^^{[0,T];T^) the set of infinitely differentiable func- 
tions in the variable x and t which are periodic in x, and we denote by C^er(T^) 



y £ "'°3 "" / \D°'u{x,y,z)\'^ dxdydz, (2.1) 
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the set of k-times differentiable periodic functions in the variable x with periodic 
domain T^. Next, we recall the Ladyzhenskaya inequality for an integrable function 
ip in two dimensions, 

ll¥'(2)IL^^(T3) < C||v.(;^)||^('(^3)||¥'W||^^f, (2.3) 

for every z £ [0, L]. Hereafter, C denotes a generic constant which may change 
from line to line. Furthermore, for all u g i/^(T^). with J,j.2 u{x, y, z) dx dy = 0, for 
every z G [0, L], we have the Poincare inequality 

||u(z)||^2(T2) < CL||u(2)|l^i(T2). (2.4) 

We recall the well-known elliptic estimate in two-dimensions, due to the Biot- 
Savart law for an incompressible vector field u, satisfying -u = and V/i xu = oj, 

\\uiz)\\^... <CpMz)\\l., (2.5) 

for every z £ [0,L] and every p e [2,oo) (see, e.g., [34] and references therein). 
Also, it is an easy exercise to show using the divergence free condition, V/i • u = 0, 
that 

Mz)\\Lf = i|V/. X u(z)||^p = V2||V„u(z)|l^p . (2.6) 

h h h 

Next, we introduce a generalization of the logarithmic inequality that was in- 
troduced in [4]. Similar inequalities were previously introduced by other authors 
in [6], [22] and [35]. These inequalities improve the Brezis-Gallouet [2] and the 
Brezis-Wainger [3] logarithmic Sobolev inequalities. We follow the same idea of the 
proof given by the authors in [4] and [19]. We give the proof in R'^, but one can 
derive the same result for periodic functions (see Corollary 2.2 below). 

Lemma 2.1. Let F : ^ R 6e a scalar valued function, F £ W^-'^^^{R^), and 
S > be given. Then, for any A £ (0,oo) and any radius R £ (0,oo), the following 
logarithmic inequality holds: 

.\F\\ 

\\F\\L-^(g2. < C^.Amax-^ l,sup. 

[ q>2 

log^ (e^ + Ci?^/(2+.)^^£Zm + ||VF||,...(^.) j . (2.7) 

where Cs,\ is an absolute positive constant that depends on 5 and A. 

Proof. Without loss of generality, we assume that ||F||^oojg2j = F{Q) ^ 0. Denote 
by Bfi the disk centered at the origin with radius R. Let </) £ C°°(R^) be a smooth 
non-negative cutoff function satisfying: 

0(0) = 1, 11011^^(^2) < 1, ||V(/.||^^(r2) < ^, suM4>)^Br. (2.8) 

Set C, = Fcf). From the solution formula of the Laplace equation in R^ we have, for 
any p>2, 

1 




2tt 



log(|y|)AC'^(y)dy| 



log(|y|)AC'^(y)dy| 



3D INVISCID PSEUDO-HASEGAWA-MIMA MODEL 



7 



By integration by parts, 



IC^(o)| = 



Jb„ y 



< 



P_ 

271 



\y\ 



+ 









Br\b, \y\ 



\cm\ < cp 



for some e E (0, R), to be chosen later in (2.11). By Holder inequality. 

y_ 

2 

y_ 

2 

L^{B„\B,) 



y| 

+ Cp 



livcil IICII 



L—^^P-^^Br^XB,) 



where a and /? satisfy ^ + ^ = -^qif, and a e (1, 2) is chosen such that 
oo. For instant, one may choose: 



|y|^ 



< 



4 + ^ (2 + ,5)(4 + <5) 



(2.9) 



This implies: 



|C^(0)| < Cspe"^ l|VCIL.+.(K.) IICII^^,'.-.,(^.) 



IIVCII 



+ Cplogi/2 

Now, we choose e £ (0, i?), small enough, such that: 



) llCir 2(2 + 5) 



= 1, 



(2.10) 



(2.11) 



thus we have. 



log 



— !|VCili2+.(^2) < R—R^ 
^log fe2 + i?V(2+5) IIVCII 



RecaU that p^/P < C for every p > 2, thus (2.10) and (2.11) yield 

i-i 

^/3(p-l)(jj2) 

log (e2 + i?^/(2+^)||VC|L2+.(^2) 



2-a -5 \ 1/P 1- 

ic(o)i<c5 i?— iicr 



(2.12) 



C 



2-a 



II ^Cll r 2 + «/„2 JICII 2(2 + a) , 



Observe that since ^ > 6 (see (2.9)) and p > 2, then /3(p - 1) > 2, /3^(^"p) < 
and {p- l)^(i-p' < p^. Then, 
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.1-^ -i?^/(^«(/3(p-l))^(i-^)' 



^3(p-l)(]g2) y^VF -JJ U^(p- l))>i?2/(/3(p-l)) 



1-J 



,9>2 Q 



ll^lli'm^) 



< max 1 1, sup | • (2.13) 

Similar argument shows that, 

\\C\\'-.L. , <C,,.i?^/(^(^+^)Vmax|l,sup^^^|. (2.14) 

Since, due to (2.9), R^R^^R^'f^ = 1 and {^f''^^^ < 1, (2.12), (2.13) and 
(2.14) win yield: 



|C(0)| < Cs., (l + logV(^^) (e^ + i?V(2+.) ||VC|1,...(».))) 



Now, we choose p > 2 such that: 



p = log (e^ + |1VC||^2+.(K.)) , (2.16) 

then we have, 

^5/(2+5) II VC||^...,«.Y^' < e. (2.17) 



IL 

Thus, wc can conclude from (2.15) and (2.17) that 

|C(0)| < CMmax|l,sup^^^^| log^ (e,^ + i?^V(2+^-) || VCIL...(^.)) . (2.18) 
Finally, notice that: 

|C(0)| = |</>(0)i^(0)| = !lF||^o„(^.), 

and that 

I|VCIIl2+''(R=^) = \\'^i<PF)\\L^+'iW'') 

< ||V0l|^^(]j,2) ||F||^2 + 5(^2j + ||VF||^2 + i(^2-| |j0||^oo(]g2) 

< I ^ + II IIl2 + ,S(r2) I . 

This proves inequality (2.7). □ 
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Corollary 2.2. Let _F : — > M 6e a scalar valued periodic function with periodic 
domain F G W^''^+'^{T'^), and S > be given. Then for any A e (0, oo) the 
following logarithmic inequality holds: 

\\F\\^^^j2) < 

m 

q>2 

where Cs,\ is an absolute positive constant that depends on 6 and X. 

Proof. Since F is a periodic function with periodic domain = [0,L]^, in the 
above proof we may choose R = 4L. Notice that for any p > 1 

IICII 

and that 



Cs,^ max I 1, snp ^^^^ log' (e^ + CLV(2+^) ||^^||^^ . (2.19) 



||VC|L.(s,,) < Uh^iB,,) \\^F\\LPiB,,) + l|V0IL.o(B,,) \\Fh.iB. 



4L I 



Thus, one can follow the steps of the proof of the previous lemma to prove inequality 
(2.19). □ 

Proposition 2.3. Lei u; e L°°(T3) n L^(W";,^'^(T3)), w 7^ 0. Let Vl = K ^ uj . Then 

h 

VhU e L°°(T3) and, 

( ||V;ia;||^=o,4„,3A 
||V;,u||^.,(T3) < C M\L-ij-) log + CL" ""^ . (2.20) 

Proof. Due to (2.5) we have, 

llV/iu(z)||^,(T^) < C(?||cj(z) 11^,(^3) <C(7L2/'?||w(z)||^c^(t3) , 

and 

l|V/iu(z)|i,^^i,4(.;^2) <C||V/iw(z)||^4(Tr3) , 

for every q G (l,oo) and every z € [0, L], where the constant C is independent of 
z. Set 



and apply the logarithmic inequality (2.19) with A = 1 and 6 — 2, we obtain 

{||cL>||^2f^3N ] / ||V/,,a;(z' 
^3/2 M I^(^)IIl-(t-^) Mog {e^ + CL^ 

for every z G [0,i]. This implies that 

(II V/iW||^oo,4™3i \ 
e' + CL' . 
II^IIl2(t3) / 



□ 
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Corollary 2.4. Let9,rje L°°{T'^)r\Lf'{Wl''^{T^)) , x u = j^iO-rj), [9,7]) ^ 0. 
Then V,,u £ L°°(T^), and 



l!V,u||^^(T-^) < - 11(0, r,)||L^(T^) log ( + ^^^^^^^-^^^^^'^ ) • (2-21) 



Proof. In this case one can choose, 



||(6',^7)||T2,T3^ 



in (2.7) and proceed as in the above proof of Proposition 2.3. □ 
Proposition 2.5. Let a; e L°°{T^), andu = K *oj. Suppose thatil; e Lf{Wl''^(T^)), 

h 

then {u ■ Vh)i' e L'^{T'^) and 

ll(u- Vh)7^|li2(T3) <CL2||w||^^(T^,)||Vft,V||i~,4(T3). (2.22) 

When uj = ^{9 - 77) with 9,r] e L°^{T^), then 

||(u.V,0^IL2(t3) <CL||(0,77)||L^(T-^)||V„^||^oe^4(T3). (2.23) 

Proof. Since u(x, y, z) = 0, by the Sobolev type estimate (2.3), the ehiptic 

regularity estimate (2.5), and the Poincare inequality (2.4) we have 

l|u(2:)lli4(T3) < C||u(z)||^2(T^) l|u(z)||^i(T3) < Ci||u(z)||^i(T3) < CL||w(z)||i2(T3) , 

for a.e. z G [0, L], Then, by Holder's inequality and the above inequality, we obtain 

||(U- Vft)V'(z)||i2(T3) < ||u(z)||^4(t3) ||Vft?/'(z)||^4(T3) , 

< CL \\U}{Z)\\12^J3^ l|VhV'(^)lli4(T3) • 

Integrating the last inequality with respect to z over [0, L] will imply: 

II (U • V/OVII L\l^) < C'^^^^ ll'^ll l2(T^) II V/iV-ll L~^*(T^) 

< C^i' ll^llL-(r') l|V,.V||i==^4(T3) . (2.24) 

When, UJ ~ jj-(9 — ri), its enough to notice that ||a;||^ooj.,jf3-) < ||(0, r/)||L=c(..j.3j. This 
concludes the proof. □ 

Corollary 2.6. Let uj e L'^{T^), andu^K *lo. Suppose that ip S Lf {W}'^ {T'^)) 

h 

and (j) e L2(T3) , then 

|((u. V,0^.,<^)| < CL^I^M\^.^^.^ \\WiMl-,\i-) U\\L\f-) (2-25) 
Proof. By Holdar's inequality and (2.24) in the proof of Proposition 2.5 we get 

|((u- Vh)V',0)| < ||(U- V/i)V'||i2(T3) 11011^2(^3) 

< II'^IIl^ct^) \\^hi^\\L^-^^ij-i) II</>I1l2(t3) . 

□ 
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3. The "Linearized" System 

Our task is to prove the global existence and uniqueness of solution of system 
(1.7). We will prove our result using an iteration procedure. For this purpose, we 
introduce and study, in this section, a "linearized" version of (1.7), about a given 
state uj, which will play a crucial role in the iteration procedure. The "linearized" 
system is considered subject to periodic boundary conditions with a basic domain 
T"^ = [0, L]"^ C M'^, over any fixed arbitrary time interval [0, T] and is given by 

^ + (u • V;,)0 - (7o^ = 0, ^ + (u.V,0r7 + ;7o^=0, (3.1a) 
ot oz at oz 

u = K (3.1b) 

h 

0(O;x) = 0°(x), 7y(0;x)=77°(x), (3.1c) 

where w is a given periodic function. 

Theorem 3.1. Let T > be given, and let Cj £ C'^{[0,T];C^^^{T^)), {9° ,ff ) € 
Cper(T^)xCper(T^)- Then (3.1) has a unique solution {9, rf) £ C2([0, T]; ^^^^(T^) x 
Cl,{T^)) with 

||(0(t),r;(t))|li^.(T3)-||(0°,7;°)||j^,(^3), (3.2) 
for all t G [0,r] and all p G [l,oo]. Moreover, 

sup ||(V„0(i),V^77(i))llL-'^r') <e''^^^^||(V,.0°,V;,77°)||j^..,4,T3., (3.3) 

sup \\{d,e{t),dMm\L\T^) < H{T)T+\\{d,9",d,if)\\^^, , , (3.4) 
*e[o,T] ^ ' 



and 
te[o,T] 



sup ||(5«0(t),5t,?(i))|!L.(T3) < H(T)r+ 11(9,00,9,7,0)11^^.(^3 



+ CLV(^)^||(V,0O^y^^O)|| 



(3.5) 



Furthermore, 



sup ||VhU(i)||^oo„3) 

te[o,T] 



< - ||(0°,.°)|L^(,3,log e- + CX-e^^(^)^ ||(.o,,0)||^.;;^ 1 • (3-6) 



Here u = * [9 — rj), J{T) and H{T) are constants that depend on the domain 

, the norms of lj and 9 as well as the time T and are specified in (3.10) and 
(3.13), respectively, while C is an absolute dimensionless positive constant. 

Proof. If (6*0, ??0) ^ 0, then the unique solution of (3.1) is {9{t),r]{t)) = 0, and 
there is nothing to prove. Therefore, we assume that {9^,rf ) ^ 0. Since uj £ 
C^([0, T]; Cpg^(T^)) then one can easily show by classical elliptic theory that u e 
C\[0,T];C^^.^{T^)). By the method of characteristics ([7], Chapter 3, [20], Chapter 
1) we can ensure the existence and the uniqueness of a classical local solution 
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(0, 7]) e C2([0, T,na.];C^,,it^) xC^^.iJ'')) of (3.1) for some short time < T^ax < 
T. Since (u, — 1) is a divergence free vector field in R"^, multiplying the evolution 
equations in (3.1a) by \0f~ and \rj\^~ rj, respectively, and integrate over T"^ yield 
that 

II^(OIIlp(T^) = ||^°ILp(T3) find |h(i)|lLP(T-') Ik/"||LP(T3) ' (^•'^) 

for any p G [1, oo) and every t € [0, Tmax]- Thus, 

ll^(i)llL-(T^) = pl™^ II^'IIlpct^) pi™, II^°ILj'(t^) II^1l°°(t3) ' (3-8a) 

h(i)llL-(T3) = pi™ hllLP(T^) = plm |h°ILp(T3) = lh°IL~(T3) ' (3-8b) 

for every t e [0,T,„a^]. The uniform estimates (3.7) and (3.8) allow us to deduce 
that the solution (9, rj) can be extended on each interval [nTmax, ('T-+l)'7maa:]n[0, T]. 
So, for every arbitrary T > 0, we have found a unique classical solution (0, r/) G 
C2([0,r];C2^^(T3) xC2g^(T3)) of (3.1) that satisfies the estimate (3.2). Since the 
solution {0{t),T]{t)) € C^eri"^^) X C'per(T^)> ^c Can take V,, of (3.1a) and multiply 
by |V/i6'pV/i6' and by |V;i77pV;i?7, respectively, and then integrate horizontally to 
obtain: 

1 dWVhdWlijT^) _ Up 6>||V/t6'Ht4(T3) 

4 dt 4 9i 

((a,u • S/h)0 ■ d^ei^hOf + {dyiL ■ S/h)0 ■ dyei^hdf) dxdy = 0, 

which yields: 

l d\\VHOit-z)\\l.^j.) t/p an V,,g(t; 2)11^4 (t3) ^ „^ 4 

4 — < ||V,u(t;z)||^oo(T3) ||V,0(t;z)||^4(T3) 

Notice that since the solution here is classical, all these estimates are rigorous and 
there is no need to justify them further. Using Proposition 2.3 together with the 
above inequality we obtain: 

"^^^'fi-'' - , ,,||,,^^^ _ ,3,, 

where 

, l|VhW(0llL~^4(T3)\ 

J(T) := C sup \\0:{t)\\^^„,. log + CL^ sup „,,^,„ . (3.10) 

te[o,T] \ te[o,T] J 

Applying the method of characteristics and then Gronwall's Lemma to (3.9), one 
can show that, 

||V,0(t; z)|li4(T3) < II V,0"(z + Upt)\\\,^^^,^ 

<e4-^^(^)^||V,0°f^^^4(^3), (3.11) 

for every t £ [0,T], and every z e [0,L]. Similar inequality holds also for 77. This 
proves estimate (3.3). Since 9 € C^([0, T]; Cpg^(T^)), we can differentiate (3.1a) 
with respect to z and get that dz9 satisfies 

+ (a,a • v,)0 + £1 • ^hdj - c/o0 = 0; (3.12) 
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and a similar equation holds also for dzTj. Since (u, —1) is a divergence free vector 
in R'^, taking the L'^iT'^) inner product of the above equation with 526' implies that, 

Therefore, by Corollary 2.6, we have 

1 d \\dzd(t)\\'\2„3s 

^ ' \te[o,T] J 

where in the last inequality we used (3.11). Similar inequality holds also for dzTj. 
Set 

7^(T) :=4Cii/2gJ(T)T||(y^^o^y^^^O)|| ||a,(I.(0IL2(T^). (3.13) 

'■'^ ^ ' te[o,T] 

Thus, 

d\\{dze,dzr])\\^2(^.^3s^ 



< H{T), (3.14) 



dt 

which implies 

ll(a.0(O,a,r;(O)IL.(T3) < II (a,0°,a,7yO) 11^.(^3) + i?(r)r, (3.15) 

for every t e [0,r]. This proves estimate (3.4). Finally, by Proposition 2.5 and 
(3.11) we get 

|l(u.V,)0(t)|l^.(T3) <CLV(^)^||V,0"||^^,4„3. sup |l^(t)||^^(T3), (3.16) 

^ > tg[o,T] 

for every t G [0,T]. Similar inequality holds also for ||(u • \7 ii)j]\\j^2f^j3y 
Next, we observe that since (3.1a) is equivalent to: 

f ^-(fi.V.). + C/„|, |?--(u-V.),-^o|, (3.17) 

then by (3.4) and (3.16), we obtain estimate (3.5). 

Set u ^ K * u), where u) = jj-{6 — tj). Hence, the logarithmic inequality (2.21) 

implies: 

C ||(V,A^),V„7y(^))||^^»^4(^r3)^ 
||V.u(t)||,.^,3,<-H(.(t),,(.))H,.^,3^1og(^e- + CL- ||(.(,),,(t))||,,;3; J 

for every t e [0,r]. Thanks to (3.2), \\{e{t),i^{t))\\^2^^3^ = || (0°, r/O) ||j^,(^3) and 

||(0(t), ?7(t))||Loc(T3) — ''°)|Il=^(t3)' every t £ [0,T], as a result (3.2) and 
(3.3) imply (3.6); and this completes the proof. □ 
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4. Global Existence and Uniqueness 

In this section, wc aim to prove global existence and uniqueness of solutions to 
system (1.7), subject to periodic boundary conditions with a basic domain T'^ = 
[0,L]'^ C M^, over any fixed arbitrary time interval [0,T]. We recall in more details 
system (1.7) 

|? + (u.V,Oe-C/o|^=0, § + (u • V,)7? + C/o^ = 0, (4.1a) 
at oz ot oz 

u = K*u, uj^^{e-r]), (4.1b) 

0(O;x) -0O(x), 7K0;x) =7?°(x). (4.1c) 
Next, we give a definition of a solution to system (4.1) or equivalently (1.7). 

Definition 4.1. Let T > 0, and {e°,ri°) e L°°(T3) with {^hO^^^h^f) e L^,f (T^), 
{d^0°,d^'q°) G L2(T3) be given. We say that (61,77) is a weak solution of (4.1), if 

(0,77) e C([0,r];L2(T3)), (V,0, V,.7;) e L°^{[0,T];LZt{T')), 

{d.e,d.v) e L°°([0,r];L2(T3)), {dte,dtv) e L°°([0,T];L2(T3)), 
61(0; x) = 6'°(x) and 7/(0; x) = 77° (x), 
such that (0, 7/) satisfy (4.1) in the weak formulation (i.e., a solution in the distribu- 
tion sense). That is, for any scalar test functions (p{t;x),i(j{t;x) £ C^j.{[0,T];T^), 
such that (/3(r;x) = 0, ■0(T;x) ^ 0, 

''^ (^(^)' ^) {ue, V,^) + Uo £ (e{s), ds = (0°(x), ^(x, 0)) , 

(^(*)' ^) - £ ^''^) - (''(^)' ^) = (^°(^)' ^(^' 0)) • 

where u ~ K * uj = 7k- K * (6 — ri). 

h h 

In the case of periodic boundary conditions, it will be sufficient to consider only 
test functions of the form 

ipit; x) = xm(i)e'"™-", m x) = 7k(^)e2"''•^ (^ ^a) 

Xm,7keC°°([o,T]), xm(r)=7k(r) = o, (4.2b) 

for m, k € , since such functions form a basis for the corresponding larger spaces 
of test functions. Thus, {6,ri) is a solution in the distribution sense of (4.1) if 

^27rim-x\ /•„\ / /,,/) y7 „27rim-x 



"'0 

f/o ^g(g), ^ Xm(s) ds ^ (^"(x), ) Xm(0), (4.3) 

(7;(s), e^-"^-) 7k(5) ds - (U77, V,e2"k-x^ ^^(,) 
Jo 

/•T / a 27rik x \ 

UoJ^ ^»7(s),^^^^^7k(s)ds = (r/"(x),e2"''-)7k(0), (4.4) 
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for all m,k G 1? . Wc notice that if (6*, 77) is a weak solution, then by Proposition 
2.5, (u • Vh)6', (u • Vh)?7 e L°°([0,r], 2.2(1^3)). Consequently, we can apply Lemma 
1.1 ([33], page 169) with X = L'^iJ'^) to conclude that: 

= ^1+ t fv„ • iue){s) + Uo^is)) ds, (4.5) 



dz 

r? = 6 + ^* (^h ■ inv){s) - Uo^Js)^ ds, (4.6) 
for some ^1,6 € ^^(T^), a.c. f G [0,r], and, 

-(t) = -V..(u0)(t) + C/o^(t), (4.7) 

^(i) = -V..(ur;)(t)-[/og(t), (4.8) 

in L2(T^), a.e. i € [0,T]. If (6*, 77) is a solution of (4.1), and equivalently (1.7), ac- 
cording to Definition 4.1, then (9^6*, 9^7?) G L°°([0,r];L2(T3)), {dte,dtr]) G L°°([0, T]; L2(T3)), 
and by Proposition 2.5, (u • V/i)e', (u • Vh)v & L°°i{0, T],L^{T^)), then we have 

^ = -V,-(u0) + ;7o|^, inL-([0,r];L2(T3)), (4.9) 

1^ = -V, • (U7,) - {/o|^, ini°°([0,r];L2(T3)). (4.IO) 
That is, for any ■0i,'i/'2 € ^^([0, T]; ^^(tS)), 

{t),^,{t)\ dt^ {-Vh ■ iuem, Mt)) dt + UoJ^ (—{t),Mt))dt, 



,^ ^g^{t),Mt)J dt^ {-Vh-{m^){t),Mt)) dt-Uo ]^ {^^^{t),Mt)J dt. 

Theorem 4.2. Let (6*°, 7?°) G C^^rClT^) x C^er{^^)> andT >0 be given. Then there 
exists a unique weak solution {9, rf) of (4.1) in the sense of Definition ^.1 such that: 

sup ||(0(t),77(t))||LP(T3) = ||(^°,^°)|Lp(t3) , (4-11) 

for any j» G [1, 00] , 

sup \\{VhO{t),Vhrim\i.-^\T^) < MT) \\{VhO\Vhif)\\j^^A„,. , (4.12) 
te[o,T] ^•'^ ^ ^ 

sup |j(a,0(t),5.77(t))||L2(T3) <i/o(T)||(a.0",a.r;")||j^. 3., (4.13) 
te[o,T] ^ ^ 

sup ||(at0(t),9ir;(O)|iL^(T3) <i?o(T)||(a,0",a,7yO)||j_. 3. 

te[o,T] ^ ' 



+ CiJo(r) ||(e°,^")||L==(T3) \\i'^he\Vu7f)\\^^A^^.^ ■ (4.14) 
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Furthermore, 

sup ||VhU(t)|l^oc„3) < 

te[o,T] 

Here u ~ * {0 — rf). C is a positive dimensionless constant, Co,Jq{T) and 

Hq(T) are constants that depend on the domain as well as the norms of the 
initial data {9'^,ri'^) specified in (4.22), (4.24) and (4.31), respectively. 

Proof. If {6'^, if) = 0, then {9{t),r]{t)) = is a solution, and by the uniqueness 
part that wc show later, it will be the only solution, and the statement of the 
theorem follows. Therefore, we assume that (0°, rp) ^ 0. The idea of the proof is 
involving the construction of a sequence of approximate solutions of (4.1) based on 
the "linearized" system, as it was introduced in section 3; and then passing to the 
limit to a weak solution. We proceed with the following steps. 

Step 1: Construction of a sequence of approximate solutions. 
We consider the sequence {(^'", ?7")}J^0' where, for n ~ 1,2,..., {9^, 711"^) is the 
unique solution of the linear system 

^ + (u"-i • V,O0 - f/o V ^ 0' ^ + (u"-i.V,)j7 + C/o^ = 0, (4.16a) 
at oz ot oz 



(961" 


-0, 


dif 


■ 

OZ 


~dt 




= K 





^- , c."-! = ^ (r-1 -,?"-!) , (4.16b) 
r(0;x) = 0°(x), 77"(0;x) = 77" (x). (4.16c) 

Since (0°, rp) S Cpg^(T^) x Cpg^(T'^), one can use induction steps to establish the ex- 
istence and the uniqueness of the sequence of solutions (0", 77") € C^([0, T]; Cpg^(T^) x 
Cpg^(T'^)) to (4.16) by virtue of Theorem 3.1. In addition, thanks to (3.2), one also 
obtains the following uniform estimate: 

||(r(t),ry"(0)|lL.(T3)-||(e°,ry°)|Lp(T3), (4-17) 
for any p E [1, 00], and all t G [0, T]. In turn, this implies that for all n = 0, 1, 2, . . ., 

sup |l.."(t)||^.(T3) < J \\{9"{t),v"{mL-iT^) - 7 ||(^°,'7°)|L=°(T3) ■ (4-18) 

Claim 4.1. For all n = 0, 1, 2, . . 

sup ||(V„r(i),V,,77"(t))||L-,4„3. < Jo(r)||(V„0°,V„r,°)||j^^,. 3,, (4.19) 

where Jo(T) is a constant that depends on L and the norms of the initial data 
{9'^,ri^), as well as the time T and is specified in (4.24), below. 

Proof of Claim ^.1: By a similar argument as in the proof of Theorem 3.1, we take 
Vh of (4.16a) and then muhiply respectively by \SI hO''?^ hd''^ and iV/iT^" pV/^ry" 
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and then integrate horizontally and obtain 

gl t^o q; <4||VhU {z)\\^^^^,^\\VhO WlUt^) 

(4.20a) 



9||V/,r?"||^4,T3^ 9||V^?7"||r4,T3^ 

(4.20b) 



By the logarithmic inequality (2.21) and the estimate (4.17), we have: 



/ ||(V,0"-l(t),V;.<-l(i))Lc..4„3) 

T \^ + ^^^^ H(.o,,o)||^.^^3, 

' (4.21) 

for some positive dimensionless constant C\ and all t € [0, T\. This implies that 



where 



^0 := f ii(^°,.°)L^(.3,iog \^^c.L- > • ^'-^'^ 

Using the above estimate, and applying the method of characteristics, followed by 
Gronwall's lemma, to (4.20), for n = 1, implies: 

for all t E [0,T]. As a result, and thanks to (4.21), we have 
l|V.u^W|L.(,3^ < ^ ||(.°,.")L.^,3^1og [e^ + C.L^e-- W^,^,.;; 

Ci ( ( ||W°.V^77°)||^^.4(T3. 

= c„(i + Y-<||(e°,^'')IL-,,,,). 

for al t € [0,T]. Again, applying the method of characteristics followed by Gron- 
wall's lemma to (4.20) for n — 2 implies 

||(V.0^(i), V.ry^(t))||,..(,3) < ||(V.e°, V.^°)||,.,4(,3) , 

for all t £ [0,T]. Next, we proceed by induction. Assume that: 

||(V,r(t), V,r;"(t))||i^^.(T3) < ||(V/.e", V,77°)||^„^4(^3) , 
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for all t e [0,T], where 



(i + D! 



Then by (4.21) we show that: 

r ( ll(Vh6'"(0,Vft77"(<))||L-,4..j,3/ 



for all t S [0;T^]- One can easily see, by the properties of the exponential and 
logarithmic functions, that 

||(V,0O,V,77O)||^^,4,^3. 
<Ci |g„(t)i + log|e2 + CiL2 -— — ^ ^ 



II(^°,'?°)IIl-(t3) 



Thus, 

|lV,,u"(t)||^,^(T-^) < 5„+i(t), 

for all t G [0,T]. Applying the method of characteristics followed by Gronwall's 
lemma to (4.20) for r? + 1 implies: 

II (V;,r+i(i),V;.77"+i (0)11^0.^4(^3) < e«"+^«* ||(V^0°, V^77")||j^^^4(^3) . (4.23) 
for aU t e [0,T]. Setting 

Jo(T) = e^»^^ ' " \ (4.24) 
thus we have proved, for all n = 0, 1, 2, . . ., 

|!(V„r(t), V„77"(i))llL-^T3) < MT) ||(V,0°, V,r7°)||L-^(T3) , (4.25) 
for all t G [0, r]. This completes the proof of the claim. 
Claim 4.2. For all n = 0, 1, 2, . . 

sup ||(9.r(0,a.77"(t))llL^(T3) < Ho{T) ||(9.0°,a.77°)||j^. 3. , (4.26) 
te[o,T] ^ ' 

where Hq{T) is a eonstant that depends on L and the norms of the initial data 
(O'^jT]'^), as well as the time T, and is specified in (4.31). 
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Proof of Claim 4-2: The evolution equations of dz9" and dztj" are, respectively, 
given by: 

An ri^fi^ 

+ (u"-i • V,Oa.r + (9,u"-i • V,0^^" - C/o^ = O, (4.27a) 

+ (u"-i • V„)a,77" + (5.u"-i • Vu)if + C/o^ = 0. (4.27b) 

When we take the L'^{T^) inner product of (4.27a) with 9^6*", and of (4.27b) with 
dzTp , we obtain: 

ld||5^6l"f,2.^^3^ , 

1 d ||9^77"||?2™3^ ^ , 

^ ^) + ((a.u"-i . VH)t,dzv'') = 0. 

Since llSza;""^ 11^2(^3^ < j^\\{dz9"~^,dzr]"~'^)\\^2^j.3y then by Corollary 2.6 and 
Claim 4.1, wc have: 



< 



dt 

■^MT) ||(V,.0", V„r;0)||^^^.(^3) ||(9.r-\a,,7"-i)||j^.(^3) , (4.28) 
for some positive dimcnsionlcss constant C2. When n = 1, we have: 

^ d^^^ ^ IT75^o(^) v.^°)||,.,(,3) ||(a.0°,a.ry°)||^2(,3) . 

After we integrate with respect to time we obtain: 
||(9.0i(t),a.r;i(i))||L.(T3) < 

II (9.0°, 9.7,0)11^^2(^3) (^l + -^Jo(r)t||(V,0O,V,77°)|L-^(T3)) ■ (4.29) 
For n = 1, the inequality (4.28) reads: 

d||(a.6'2,a.7?2)|| (J „ n n „ „ 1 1 „ 

dT^^ - 1^75^°^^) II(^'«^°'^'^^°)IIl-(t3) ||(5.^\a.,yi)|L.(,3) . 

Thus, using (4.29) and integrating with respect to time, we reach 



^ 1 f C2 ^'^ 



||(9.^^a.,0)||^2(,3)Ef[(fT7J'^0(^)^IK^'«^^^'«^°)|lL-(T 

fe=o ■ ^ 

Wc proceed by induction. Suppose that: 

||(5,0^a.77O)||^^2(T3)E fci TT72-^o(^)^I|(^'^^"'^'^''°)|Il--(^^^ ■ (4-30) 
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After we substitute (4.30) in inequality (4.28), for n+ 1, and integrate with respect 
to time wc have, 

||(9.r+i(t),9.,y"+i(0)|L.„3,< 



r; 

|(a.0^a.ryO)||^.(,3)E^(f^72'^o(^)^||(V.0^V.r;°)||,_^,(,^ 



ri+1 -I / \ fc 



fc=0 

Set, 

then we have just proven that for n = 0, 1, 2, . . ., 

||(a.0"+i(<),a,7f+i(t))||j^.(^3) < Ho(T) ||(a.0°,a.7y°)||L.(^3) , (4.32) 

for all t £ [0,T]. This concludes the proof of Claim 4.2. Finally, by Proposition 2.5 
and estimate (4.17), for all n = 0, 1, 2, . . ., 



||(u"-^ • V,)r(t)||^.(^3) + ||(u"-^ • V,)rr(t)|L2(T3) 

< Ci||(0°,r;")||j^^(^3) ||(V;,r(t),V,r;"(t))||i^oo^,4(T3). (4.33) 
As (6'",77") is a classical solution then by (4.26) and (4.33) we have 

sup \\{dte-^{tidtifmkHt^)<MT)\\{d.e\d.rf)\\^. , 

te[o,T] ^ ' 

+ CLMT)\\{e\if)\\^^^^,^ \\{'^he\Vhv'')\\^^Afj^y (4.34) 

Step 2: Convergence of the sequence {(6*", ?y")}^o- 

Wc will show that the sequence {(0", ?7")}5^o a Cauchy sequence in C([0, T]; L^(T'^)). 
Claim 4.3. For every n = 0, 1, 2, . . . , 

\\(an+l _ nn „n+l _ n\|| ^ 

II a ,ri r] ^11^2(^3^ ^ 

1 / 2C 

IK^^ IIl^(t3) (lT7l^o(r)^ IK^'^^"' V.^«) 11,^,(^3) 

/or some positive dimensionless constant C3. 

Proo/ 0/ Claim 4.3: Wc notice that {9''+^ - 9'') satisfies 

^i^:^^^ + (u« . v.)r+i - (u"-i . v,o^" - = 0, 

u"-^ = K * u}"-\ u" = i^: * w", a;"-i = —{9"-^ - 77"-^), = — (r - r;"), 
?i h 2L 2L 

r+HO;x) = 6'"(0;x) = 6l°(x). 

This can be rewritten as: 

^"^ + ((u" - u"-i) . v,)r+i + (u"-i . v,)(0"+i - 0") - f/„^(^::^_^ ^ 0, 

and, taking the L^(T^) inner product with (0"+^ — 0") yields: 



ld||0" + l -0"||%„3, 



2 
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Thus, by Corollary 2.6 we have, 



d 11^^^^^ 0^^ II ^ 

< |(((u"-u"-i)-V,06'"+\6i"+i-6i")| 



2 dt 

< CgL^/^ II V^^'^+l 11^00^4(^3) ||tj" - w""^||^2(t3) ||6'""^^-^'"IL2(t3) 

for some positive dimensionless constant C3. Similar argument with — ry" will 
show, 

1 d||?7"+i-r;"||%,„3, 

" ' ll-f- (T ) ^ ^ rl/2 IIy7 /in+lll II n n-l\\ \\ n+1 n\\ 

2 Jt <<~^3^' W^ht^ ||L~-^(T^)lr IIl = (T^)II^ IIl2(T^) 

Il2(t3) - xlK""""" %??"-?7" 'J||l2(T' 
(4.19), we have: 

a|Hf (7 /IIl^(t^) 
Jt 

< -^if„(r)||(V,0°,V,r;«)||j^^^.(^3jc."(<)-c."-i(O|L.(T3) 

^ :^^^'"(^) IK^''^"'^''^''')IIl-^(t3) IK^" - - '/"~')|Il=(t3) ■ (4-36) 

Notice that 6'"+i(0;x) - 0"(O;x) = 77"+i(0;x) - r/"(0;x) = 0, and notice that 
11(6*1 -6*°, -7;0)||j^2(t3) < 2 11(610,770)11^^2(^3) V71 = 0,1,2,.... Integrating inequal- 
ity (4.36), for 71 = 1, with respect to time implies: 

II i9'it) - 6\t),Tf{t) - r?i(0) L.(T3) < JT^Ko{T)t \\ (V,.0O, V,.7;") ||j^^^.(^3) || (^", ^?") |L^(t3) 
for all t e [0, T]. Estimate (4.36), for n = 2, together with the above inequality give 

rf||(^3_^2^^3_^2)||^^^^^^^ 



< 



dt 

-^ifo(r) ||(v^0°,v,,77O)||j^^^.(^3) ||(^' -??^)||l^(t3) 

2 



< II(^^'?°)IIl^(t3) (|§^^o(t)||(v,0",v,7,°)||,^,(,3)) t. 

Integrating the inequality yields: 

11(03 _^2^^3_, ^2)11^^^^^^^ < i||(e",^°)L.(,3) (^|g|A'o(r)i||(V.0",V.7;")||,.,(,3))' 



Proceeding by induction, suppose that 
||(6i"+i-r,7?"+i-77")' 



< IK^"''^")|L.(T3) ( ;^i^o(r)t ||(V,e",V,77'^)||,^,(,3) 



(4.37) 
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Then by (4.36), and the induction assumption (4.37) we obtain 

Jt 

< CMT) ||(V,^°,V,77")||j^^^.(^3) ||(0"+i ~,7")||^.^^3) 

< i II (^Ko{T) II (V.0°, V.^°) |L .,,,3,) r. 



n 



Integrating the inequahty imphes: 

ii(0"+'-^"+\'7"+'-^"+')IL^(t3) 

^^^.W^^'^'^'^IW^) (|§^o(r)t||(v.0°,v.^°)||,.,.(,3))"^'. 

This, in turn, proves Claim 4.3. Now, let m > n, we conclude from Claim 4.3 that: 

sup \\(e"\t)~e-{t),rr{t)^rrm\^. . 

te[o,T] 

<I|(^°'^°)|Il^(t3) E ^(|§^o(T)T||(V.0°,V.^O)||,.,(,3 

k=n+l ^ 

This implies that {(fi'",??")}^^! is a Cauchy sequence in C{[0,T];]L'^{T^)) since 

fc=o ■ ^ / 
Then there exists (6',77) e C{[0,T];L,'^{T^)) such that: 

(r,7n^(^,'7) in C([0,T];L2(T3)), (4.38) 
as 71 —> 00; and hence, 

||(0(i),r/(t)||L^(T3) = hm |l(r(i),r?"(t))|lj^.(T3) = \\{0'^ (4.39) 

for all t E [0, T]. Define u -.^ K * co = jj^K * (9 — rj), then by the elliptic regularity 
estimate (2.5), 

||u"(t; z) - u(i; z)||i2(T3) < CL ||c."(i; ^) - uj{t; z)\\l2^^s^ 

< C |l(r (t; z) - e{t- z), 77"(t; z) - ^{t- z))\\l.^^^.^ , 
for every n = 1, 2, ... After we integrate with respect to z over [0, L] we obtain: 

||u"(t) - u(t)||^.(T3) <C\Wit) - 6{t),,f{t) - 77(t))|lL^(T3) , (4.40) 
for all 71 = 1, 2, . . .. Thus, we conclude that 

u"^u in C([0,r];L2(T3)), (4.41) 

as 71 — >• 00, and by (4.40) 

sup ||u(<)IL2(t3) < C sup ||(0(t),7y(t))||i^.(T3)=C||(0",7?O)||j^. 3.. (4.42) 

Step 3: Passing in the limit to a weak solution {0,r}). 

For each n = 0,1,2,..., {9^\ri") is a classical solution of (4.16), thus (^",77") is 
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also a weak solution of (4.16) in the sense of Definition 4.1, then for any Xm,7k G 
C°°([0,T]), with XmiT) - and 7k(r) = 0, 



Jo 

/■^ / a 27rim x\ 

+ [/o (e"{s), ^ ds = ) Xm(0), 



Jo 



/•J / a 27rikx\ 

t/o 7k(5) d5 = (r;", e^"''-) 7k(0). 



for all m,k e Z^. The strong convergence (4.38), obtained in Step 2 of the proof, 
implies that 

(r(s),e2— )x;.(s)ds^ r {9{s),e'^^^--)x^{s)ds, (4.43) 



(77"(s),e2-k-)7k(s)ds^ / (rKs),e2-''-)7k(s)ds, (4.44) 



^"(g), JXm{s)ds^ y^'^ dz /^'"(")^"' (4-4^) 

/ a 27rik x \ /"^ / «p27rik x \ 

^V'(s),^^^7k(s)ds^ i, (^/W,^^)7k(s)ds, (4.46) 

as n — > oo. It remains to show the convergence for the nonlinear terms. By adding 
and subtracting (u"^^6', V^e^'^™ '') and (u"^^?^, V/iC^'^*'^ ''), respectively, we get 

r (u"-ir , V^e^-™-) xm(s) ds - r {u9, V,e2™) Xm{s) ds = 
Jo Jo 

r (u"-i(r - 0), V„e2™) xm{s) ds + r ((u"-i - u)r , V^e2™) Xm(s) rfs, 
Jo Jo 

and, 

V7 „27rik-x\ j„ / /,,^ V7 „27r2k 







(u"-i77", V^.e^-'^-) 7k(s) - / (u77, V/.e^"''-) 7k(s) ds = 

Jo 

(u"-i(7?"-r,),V^e2-*-)7k(s)ds+ / ((u"-i - u)r , V/.e^-^'^-'') 7k(s) ds. 



By Holder's inequality, we have: 

i-T 



f (u"-i(r-0),V^e2™)xm(s)ds 
Jo 



< 



I-T 



V,e^— 1|^^(^3) ||u"-^(s)||^.(^3) \\e-it) - em^.^^s^ |x™(s)| 



ds. 
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Since |xm(i)l < C, for every t e [0,r], and || Vhe^^^-^H^^^^^a) < f - t^^en from the 
above and (4.42) we obtain 



T 

27rim-x\ 



V')|L^(t3) sup w\t)-e{t)\\ 



L " "" ' te[o,T] 



Similar argument, and the use of (4.39), wiU imply 

C 



Jo 

We can conclude, similarly, that 

(u"-i(7y"-?/),V;,e2'^*'^-^)7k(.s)ds 



<-T\\{e'\rf)\y sup ||u"-^(i)-u(t) 
^ ^ ' te[o,T] 





and 

T 



<T^I|('^"'^")|Il^(t3) sup 

^ ^ ' tG[0,T] 



Jo 



<^r||(0O,ryO)||j^. 3, sup ||u"-i (t) - u(t) | 
^ > te[o,T] 



Thus the convergence, 

(u"'i(r - 0), V^e^"-^) xm(s) ds ^ 0, (4.47) 







r ((u"-i - u)r , V,e2™) xm(s) ds 0, (4.48) 
Jo 



T 



{u"-\rf - 7?),V;,e2"'"'')7k(s)ds 0, (4.49) 



r {{n'"-' - u)77", V/.e^"''-) 7k(s) ds ^ 0, 
Jo 



(4.50) 



as n oo, follows immediately by (4.38) and (4.41). The convergence in (4.43), 
(4.45), (4.44), (4.46), (4.47), (4.48), (4.49) and (4.50) show that {0,r)) satisfies: 

Jo Jo 

+ C/oy^ {B{s), J xm(s) ds = (g", e^""*-) xm(0), (4.51) 



T 



(rKs), e^-k-x^ ^^(^) _ / ^^^^ y^e^^i^-^) j^{s) ds 
Jo 

fT j Q 27r2k-x ^ 



-f/oy^ ^r?(s),^— ^7k(s)ds = (7yO,e2-'^-'')7k(0), (4.52) 
where u = * w = 57 A' * (^^ — r/). 
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Now, we will check that 9{0) ~ 9^ and 7^(0) = ■tf. As we showed before the 
statement of Theorem 4.2, for any V'i,V'2 e L^{[Q,T]; L'^{T^)), 







do 



dt 



(-V;,-(u0)(t),Vi(t)) dt + UoJ^ (^^^{t),Mt)) dt, 



drj 

dz 



{t),Mt)) dt 



Choosing ipi{t; x) 



„27rzk-x 



7k (i) and using integration by 



parts and integrating with respect to time imply that 



') Xn.{s) ds - r (ue, V.e^"--) xmis) ds 
Jo 

Xm(s)ds = (^(0), 



dz 



Ox. 



and 

fT 



(vis), 



) 7k(s) ds - r (urj, V„e2-k-x) ^^(,) ds 
Jo 



C/n 



77(s), 



de 



27rik-x 



dz 



7k(s) ds = (?7(0), 



^27r2k-x 



)7k(0). 



By comparison with (4.51) and (4.52), we see that 



)Xm(0)-(^^^e2™)Xm(0), 



(77(0),e2-''-)7k(0) - (77",e2^^''--)7k(0). 

for all m, k G and any function Xm and 7k. We can choose Xmj7k such that 
Xm(0),7k(0) 7^ and therefore 

(61(0) - 6*", e^^™-^) = (7y(0) - e^^^-x^ ^ q 

This implies 

6'(0) = 6'°, ?7(0) = 77°. (4.53) 

By (4.17), the sequence {(6*", ry")}^^^ is uniformly bounded in L°°([0, T];'Lp{T^)). 
Estimate (4.19) shows that {(V^^", V^r;")}^^o is uniformly bounded in L°°([0, T]; 
L^^''(T3)). Estimate (4.26) shows that {(9^6*", a^?7")}^=o is uniformly bounded in 
L°°([0,r];L2(T3)), and the estimate (4.34) shows that {{dte"- , dtif} is uniformly 
bounded in L°°([0, T]; L^(T'^)). We can pass, by the Banach-Alaoglu theorem, 
to a weak-* convergent subsequence (that we will still denote by {{9" , v")}'i^=o) i'^ 
i°°([0,r];L°°(T3))such that {(V/.^, V,,7?")}^^o is weak-* in i°°([0, T]; L^^^(T3)), 
{((9^e'",a^77")}^^(, and {{dtO'' , dtT]'')}^^^ are also weak-* convergent in L°°{[0,T]; 
L^(T^)). Then, by the uniqueness of the limit, we have: 



9", 77") 



^ (0,77) in L-([0,r] 

(V^r , V„77") - {WhO, WhV) in i°°([0, T] 

(a.r,a,77") - (^.0,9.77) in L°°([o,r] 

(at0",at7,")-(at0,at77) in L°°([o,r] 



Lf(T^)), 

L2(T3)), 
L2(T3)), 
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as n — )■ cx) for any p G [1, oo]. Then, it follows that: 

sup ||(^W,'7W)I1l-(T^) < 
te[o,T] ^ ' 

for any p S [1, oo], and by (4.53) we conclude that 

sup \\{e{t),r]{t))\\^p„3. = I|(^°''?*')|Il''(t3) ' 
te[o,T] ^ ' 

for any p G [1, oo]. Also, 

sup ||(V„0(i), V^?7(i))||L~.4(T3) < Jo(T) ||(V^e°, V^r/°)||j^oo,4 3. , 
te[o,T] ^ ^ 

sup ||(9,^(i),a,r;(0)|lL=(T3) < Ho{T) \\{d,e°,d,v'')\\^, , , 
te[o,T] ^ ' 

sup 11(3*0(0, atr;(i))||L2(T3) < i/o(T) 11(9,0", a,77°)L^„3) 
te[o,T] ^ ^ 

+ CLJo(T) ||(0",r;O)||j^^(^3) \\i^h0°,^hv')L^A^j.~, ■ 

We just showed that {6,ri) is a weak solution of (4.1) in the sense of Definition 
4.1. Finally, by a similar argument as in the proofs of Theorem 3.1 and Proposition 
2.3, we have 

C I ||(V„0O,V„7yO)||j^^,. 3) 
sup ||V,u(i)||^.(T3) < - (0O,ry°) 3)log e' + CL^MT) 

The uniqueness of the weak solution (9, -q) is an immediate consequence of the next 
Corollary 4.3. This completes the proof. □ 

Corollary 4.3. Let (0?, 77?), (^g, %°) G ^^(TS), (V;,^?, V.t??), (V;,0O, V,7?0) G K,h(^'') 
and assume the existence of corresponding two weak solutions (^ijjyi), (02,'72) of 
(4.1), such that the weak solution (02, ^72) satisfies 

sup ||(Vft02(O, V,,772(0)|iL~.^T^) < J0,2 ||(V/,0^, V,,7?°)|| „o,4™3^ . (4.54) 
t6[0,T] 2,h ^ ' 

for some positive constant Jo. 2- Then, 

sup ||(0i(t)-02(t),mW-^?2(t))|lL2(T3) 
tG[0,T] 

Proof. Let us assume that (0i,r]i) and (02,??2) are two weak solutions with initial 
data, (0^,77°) and (02, '72)7 respectively, then the differences, = 0i — 02 and r/ = 
?/i — ?/2 satisfy 

— + V,, • (ui0i) - V,, • (U202) - C/o^ = 0, (4.56a) 
^ + V,, • (ui?7i) - V,, • (U2772) + Uo^ = 0, (4.56b) 
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in L°°([0,T];i2(T3))^ where ui = -^K * (6*1 - r/i) and Ua * (6*2 - 772)- 

Recall that O.r] e C([0, T]; i^crS)) and |f , ^ € L°°([0, T]; 2.2(^3))^ Thanks to 
Lemma 1.2 ([33], page 176), the equations below hold 

We can take a.e. in t the ^^(T^) inner product of (4.56) with 9{t) and ri{t), 
respectively, and get 

^ ' J^"^ =-((u.V,)g2(t),g(t)) 

< CL^'^ ||V„02(i)|li^,,.(T3) lh(OllL^(T^) ll^(i)llL^(T3) , (4.58) 

-{{{n{t))-VMtUm 

< CL^/^ W^hmimL^-^^f^) II^(^)IIl^(t3) MMl^t^) > (4-59) 
a.e. t £ [0, T] where u = Ui — U2 and w = oji — ^2- Recall that by the use of (2.5), 

||w(t)IL2(T3) < i||u(t)||^.(T3) 11(^(0, ^(0)I1l^(t3), 

for any t G [0,r]. Therefore, summing the inequalities (4.58) and (4.59) and using 
the assmuption (4.54) yield 



dt 



d(\\oml^ij^) + Mt)\\ 



2 



dt 

C 

a.e. t g [0, T\. Integrating the inequality with respect to time and using Gronwall's 
lemma, and the fact that norms are equivalent in finite dimensional spaces, we have, 

||(^l(i)-^2(i),^l(i)-^2(0)llL^(T^) 



.3-, , 



a.e. is [0,r]. This proves the corollary. □ 

Theorem 4.4. Let T > 0, let [9^ .jf) G L°°(T3) with (VhO^^Vhif) e L^,f (T^), 
(dzd'^ ,821^) € l-F'iT'^) be given. Then system (4.1) has a unique weak solution {9, if) 
in the sense of Definition that satisfies estimates (4.11)-(4.14) Furthermore, 

sup ||V,^U(<)||^=,.J.3^ < 

tfE[0,T] 

C ( ||(V„0O,V„77")||j^.„,. 3, 

-||(.°,,°)|U,,3,log e2 + CL2jo(T) iK^o^^O)!!^,;; I ' ^'-''^ 




where u = A' * a; with uj = ^{9 — rj) . C is a constant that depends only on the 

domain , Jq (T) , depends on T and on the domain and the norms of the initial 
data{9°,rp) and are specified in (4.24). Moreover, if hSi,VhVi), i"^ h02,y hv") & 



28 CHONGSHENG GAG, ASEEL FARHAT, AND EDRISS S. TITI 

^Thi"^^)' i^zOijdzrii), {dz92,dzr]2) E L^(T'''), then the corresponding weak solu- 
tions (^1, (^2, ^72) of system (4.1) satisfy: 

sup ||(0i(t)-02(t),??l(i)-?72W)|lL^(T3) (4.61) 

te[o,T] 

Proof. The proof is based again on the idea of construction of a sequence of approx- 
imate weak solutions {(^^, T?^)}ee(o,L/4] of (4-1) and then passing to the Umit to a 
weak solution of (4.1). Without loss of generality, we will assume that (0°, 77°) 7^ 0. 

Step 1: Construction of a sequence of approximate solutions {{9'^, 77'^)}jg(o ij- 
Let p{x,y,z) = ({x,y)£,{z) be a C°°-compactly supported mollifier, such that 
C(2^j y) ^ a-iid ^(2:) > are C°°-compactly supported moUificrs, with ^(z) dz = 
1 and J^^ Cix,y)dxdy = 1. Let p^(x,y,a;) = C{x,y)^''{z), where C{x,y) = 
^C(f, f ) and r(^) = )• By Theorem 6 ([7], page 630), (0°(x)*p^ r,0(x)*p-') e 
C;^^(T3) X C^^{T^)\fe € (0, |]. By Young's convolution inequality: 

\\0° * P^||i=e(T[.3) < ||6'°||iOO(T3) : ||^° * /'^||l'=°(T3) - lh°|lL~(T3) ' 



I * pIL-^t-^) = II (v/.e° * c) (^)IL4(T3) * r 



, , < ||V/i6l°|L^,4™3N , 



01 



4 /rTp3 \ 



By Theorem 4.2, and the above estimates, there exists a weak solution (0^,77^) of 
the system 

^ + (u- • V;,)r - C/o^ = 0, i^ + (u-.v^)77^ + C/o^-0, (4.63a) 

= A' * w% = -5-(r - r;=), (4.63b) 
/i 2/7 

r (0; x) = 61° (x) * 6l"•^ ?f (0; x) 77" (x) * p'^ =: r;°'^ (4.63c) 

h h 

which enjoy the foUowing estimates, 

sup \W{tW mWi^^T-) - ||(^°,^?°)|Il-(t3) . (4-64a) 

for any p S [1, 00], 

sup \\{yHe'{t),^hV'mk--\T-) < MT) ||(V,0O,V,77O)||j^^,4 3. , (4.64b) 
te[o,T] ''''' ^ ' 

sup ||(a.r(t),9.r?^(i))lk.(T3) < ffo(T) ||(a.0°,a.77")||L2„3) , (4.64c) 

tG[0,T] ^ ^ 

sup ||(atr(t),atrf (t))||L.(T3) < i?o(r) ||(a.0°,a,r7O)||j^. 3. 

*e[o,T] ^ ' 

+ CLJo(r) II (0°, ^7°) ||l==(t3) ||(Vhe°, V,77°) ||j^^^4(^3) , (4.64d) 
where Jo{T), Ho{T) arc constants specified in (4.24) and (4.31), respectively. 
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Step 2: The convergence of the sequence {(^^, ??^)}eg(o,i]j and passing in the 
hmit to a weak solution (6*, 7y) as e — > 0. 

Claim 4.4. T/ie sequence {{6" ,r]'')}^^^f) L] is Cauchy m C{[0,T];L^{T^)). 

Proof of Claim 44: Let j > 62 > Si > 0, and let {0"^ ,rf^),{e''\r]''^) be the 
corresponding weak solutions of system (4.63). By Corollary 4.3 we have: 

sup ue^Ht) - o'Ht),v''{t) - v'^myit^) < 

te[o,T] 

where Jo{T) is specified in (4.24). Since (6*°, 7/°) e iJ(T^), then by Theorem 6 ([7], 
page 630): 

(0"'^^r,°'') ^ {e",r,°) inL2(T3), 
as e ^ 0. Thus, as ei, £2 — , one can see that 

sup ||(r^(t)-rHf),r?-HO-r'=W)llL^(r^)^0, (4.65) 

t£[0,T] 

which implies that the sequence {(0^, 77^)}£g(o,i] is Cauchy in C([0, T]; L^(T'^)). 
Thus there exists (6*, 77) S C{[0,T];L,'^{T^)) such that: 

(^,77^)^(0,77) inC([0,r];L2(T3)), (4.66) 

as £ 0. Define uj = jji^{6 — rj), u = K * uj, then: 

u^^u mC{[0,T];L^{T^)), (4.67) 

as £ — > 0. Following the argument in Step 3 of the proof of Theorem 4.2, one 
can show that (^,77) satisfies (4.1) in the weak formulation. The uniform bounds 
(4.64) allow us to use the Banach Alaoglu compactness theorem and extract a 
weak-* convergent subsequence, as in the proof of Theorem 4.2, that inherits all 
the uniform bounds (4.64). And so, {9,1]) is a weak solution of (4.1) and by the 
uniqueness of the limit, the weak-* will be again (^,77). Moreover, by Proposition 
2.3 and by the argument introduced in Theorem 4.2, 

sup ||V^u(t)||^^„^) < 
te[o,T] 

C ( ||(V/.0O,V,77O)||j^^,4 3/ 

step 3: Continuous dependence on the initial data and uniqueness of the solution. 
By Corollary 4.3, 

sup \\{ei{t) - e2{t),T]l{t) - T]2it))\\j^2„3) 

te[o,T] 

llk-^o(T)T|l(V..«,V..«)|| , 
< Oe ||(,f^i - ^^2:'?! - %)||l2(t3) ■ V4-D») 

Uniqueness is an immediate consequence of the above inequality. This completes 
the proof. □ 
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Now, we consider the system: 

w + 77), u = ^{e-r,), (4.69b) 

V,i X u = uj, Vh. • u = 0. (4.69c) 

Notice that for any space norm we have: 

\\iw,LLj)\\ < \\{e,f])\\<2\\{w,LLj)\\. (4.70) 

The next Theorem foUows immediately from Theorem 4.4 and (4.70), 

Theorem 4.5. Let {w°,uj°) G L°°{T^), {VhW° ,VhUJ°) € L^A^(T3), {d^w°,d,uj°) S 
L2(T3), then the system (4.69) has a unique weak solution {w, w) in a sense similar 
to Definition 4-1, such that: 

sup ||(u;(t),La.(0)||L-(T^) <2||(«;°,ia.0)||j^^ 3^, (4.71) 
te[o,T] ^ ' 

sup ||(V,,u;(t),iV,,t^(t))||L=.4„3) < 2Jo(T) ||(Vfti(;0,LV,,cjO)|| 4.^,3^ , (4.72) 

ie[0,T] ' 

sup ||(a,u.(i),L9,w(t))||L2(T3) <2Ho(r)||(9,«;0,La,a;°)|L.„3., (4.73) 
te[o.T] ^ ' 

sup ||(atu;(i),i5tCj(0)|iL2„3) < 2i/o(r) ||(a.u;°,La.w°)|| . 3^ 
te[o,T] ^ ' 

+ 4CL Jo(r) II («;", Lc^O) IIj^^ ^^3^ || (V,,u;°, LV,,w°) ||j^o„^4(^3) ■ (4.74) 

Furthermore, 

sup ||V;,u(t)||^oo(T3) < 

*e[o,T] 



4lK-°'^-°)|lL-(T3)log I e- + 2CL-Jo(r) - |,„o. ,^o,|ur ^' ^ I ' ^''''^ 



z,h ^ ' 

where w = K ^ uj . Cq is a constant that depends on the domain and the 

h 

norms of the initial data (u''',a;*') specified in (4.22). C is a positive dimen- 
sionless constant , Jq{T)^Hq{T) depends on T and the on domain T'^ and the 
norms of the initial data {uP,uj'^) specified in (4.24) and (4.31), respectively. More- 
over, if {wlu:l),{wlLol) G L-(T3), (V,«;?, V;,.;?), {^uwI^hujI) e L^'(T3), 
{dzWi, dzUJi), {dzW2, dzU!2) € L^(T'^), then the corresponding weak solutions {wi, uji), 
(^2,^2) of system (4.1) satisfy: 

sup \\{wi{t)-W2{t),L{ui{t)-UJ2{t)))\\j^2,j3) < 

te[o,T] 

2Ce -."^ ' ||(wi - 'u;2,i(wi - W2))||l2(t3) ■ (4-76) 
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